Non-invasive characterization of transverse beam emittance of electrons from a 
laser-plasma wakefield accelerator in the bubble regime using betatron x-ray radiation 
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We propose and use a technique to measure the transverse emittance of a laser-wakefield accel- 
erated beam of relativistic electrons. The technique is based on the simultaneous measurements 
of the electron beam divergence given by v±/v\\, the measured longitudinal spectrum 7|| and the 
transverse electron bunch size in the bubble r± . The latter is obtained via the measurement of the 
source size of the x-rays emitted by the accelerating electron bunch in the bubble. We measure a 
normalised RMS beam transverse emittance < 0.5 tt mm mrad as an upper limit for a spatially 
gaussian, spectrally quasi-monoenergetic electron beam with 230 MeV energy in agreement with 
numerical modeling and analytic theory in the bubble regime. 



PACS numbers: 52.38.-r, 52.38.Kd, 52.25.Jm 

A laser wakefield accelerator [T] uses an intense laser 
pulse to generate a plasma wave, known as a wakefield, 
with a phase velocity close to the speed of light. If a 
relativistic electron beam can be inserted co-propagating 
with the wakefield, it can remain in an accelerating phase 
of the wave for a long distance, and extract substantial 
energy from the longitudinal electric field. Since plasma 
waves can support electric fields many orders of magni- 
tude stronger than those in radio-frequency (RF) cavi- 
ties, plasma accelerators provide a compact alternative 
to conventional RF linear accelerators in the future. 

In the nonlinear regime, the laser pulse expels electrons 
from its focal volume [2j but has a negligible effect on the 
remaining ions, resulting in a bare ion cavity with strong 
electromagnetic fields. Electrons from the periphery of 
the cavity are pulled inwards by the strong Coulomb at- 
traction and form a thin high-density sheath around an 
approximately spherical "bubble" O [4] . This has ideal 
linear accelerating and focusing properties for electrons 
within the cavity and at the point where the sheaths cross 
at the rear of the bubble, the field is particularly strong. 
Here, electrons from the plasma can be accelerated to 
the phase velocity of the bubble in a time shorter than 
their crossing time and are therefore trapped. Under cer- 
tain conditions, the self-trapped electron beam can have 
a quasi-monoenergetic energy distribution |5H7j. 

Most applications require the electron beam to be 
transported over long distances and/or to be focused into 
a small space with a minimum of divergence, for exam- 
ple to form high-resolution images. A good measure of 
the quality of a particle beam is its density in 6D phase 
space {x^y^ z^Px^Py^Pz)^ which is a characteristic of the 
beam. The Liouville theorem states that for noninter- 



acting particles, subject to a Hamiltonian (the electro- 
magnetic fields of the plasma) the extent of the beam in 
phase space T = J pdr^ called emittance, is conserved. 
For practical purposes the emittance is split into 2D sub- 
spaces x^Px^ y^Py and z^Pz. If the area occupied by the 
beam in the {r±^p±) plane is T^ with ±= x, then the 
normalised or invariant (transverse) emittance is defined 

:= r^/(7rmec) (1) 

where rrie is the electron mass and c is the speed of light. 
The gradients of the trajectory in the {r±^z) plane are 
more easily measured. This motivates the introduction 
of the trace space (r^,r^) where = dr±/dz — v±/v\\. 
If the area occupied by the beam in trace space is A^, 
then the geometric emittance is defined as 

:= A^/n (2) 

One can show that normalised and geometric emittance 
are related via e:^ = e^^Pz = iPz J r^j_dr± < jf3zr'j_r± 
where ^z = ^z/c- The measurement of the normalised 
emittance reduces to a measurement of the beam en- 
ergy 7, the transverse extend r_\_ and divergence of the 
beam r^, where /^^ c:^ 1 is assumed for relativistic beams. 

The extent in -space of a beam is relatively easy to 
obtain in general, e.g. by imaging its size with a scintil- 
lating screen. The extent in gradient r^-space of a beam 
is more difficult to determine in general. For a lower en- 
ergy electron beam, one standard measurement technique 
is the "pepper-pot" technique [H [9] . A lead mask con- 
sisting of a series of small apertures, significantly smaller 
than the incident beam profile, allows the extent in r±- 
space to be measured and to be related to the angular 
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divergence of the resulting beamlets at the aperture 
by assuming bahistic trajectories. This technique was 
apphed to 100 MeV electron beams from a laser wake- 
field accelerator p!QHT2] measuring normalised transverse 
emittances of order tt mm mrad. However it is difficult 
to extend this method to the GeV energy monoenergetic 
beams of recent laser wakefield experiments [131 [14] ^is, in 
addition to the difficulty of stopping GeV electrons, the 
spatial profile of the beam is small and therefore difficult 
to aperture. Another disadvantage of the "pepper-pot" 
technique is the difficulty of measuring the electron en- 
ergy 7 simultaneously with (r^, r^). The precision of the 
"pepper-pot" technique is ultimately limited by conflict- 
ing design considerations [15 . 

In this Letter, we propose and use a technique that 
allows the simultaneous measurement of (7,rx,r^) and 
hence determination of an upper limit to the normalised 
emittance on a single shot. The electron energy 7 
is measured with a spectrometer, the divergence r'j_ is 
obtained from the distribution in the direction perpen- 
dicular to the dispersion and the beam size r± is inferred 
from the betatron x-ray source size emitted by the elec- 
tron bunch. Using this technique, we show that the nor- 
malised emittance of a 230 MeV electron beam from a 
LWFA can be < 0.5 tt mm mrad and agrees with the den- 
sity scaling in the bubble regime. Our method will work 
without modification to measure the emittance of multi- 
GeV beams anticipated in the near future. 

The experiments were performed at the Hercules 
laser of the Center for Ultrafast Optical Science at the 
University of Michigan, Ann Arbor. Horizontally, lin- 
early polarised pulses with a central wavelength of Al = 
800 nm, a Gaussian full- width half maximum (FWHM) 
pulse duration of II = (32 ± 2) fs and an energy of 
Sl = (2.2 ± 0.1) J were focused to a peak intensity of 
lo = (2.0 ± 0.4) X 10^^ Wcm-^ or an ao = (3.0 ± 0.3) 
with an off-axis parabolic mirror of focal length / = 1 m 
and an F- number of F = 10. The transverse profile of 
the laser in vacuum yielded a FWHM focal diameter of 
dpwHM = (10.8 ± 0.5) /im. Typically, (27 ± 2) % of the 
laser energy is within the diameter dpwHM' Helium gas 
was pulsed through a 3 mm diameter nozzle creating a 
supersonic jet. The density profile measured by inter- 
ferometry has a FWHM length of (2.7 ± 0.2) mm and 
consists of a central plateau, which is constant to within 
10%, and linear gradients at the edges. 

Figure [1^ shows a schematic of the experimental setup. 
A permanent magnet deflects the electrons from the laser 
axis and disperses them horizontally according to their 
energy onto a scintillator screen, which is imaged with a 
CCD camera to obtain the electron spectrum and elec- 
tron beam divergence. Electrons with energies between 
50 and 500 MeV and divergence angles between ±10 mrad 
can be detected. Electron beams are observed for elec- 
tron densities between 3 and 8 x lO^^cm"^. A sample 
electron beam is shown in figure [2^. The average beam 




FIG. 1: Schematic of the (a) experimental setup and (b) inside 
of the plasma wakefield. 



charge is found to be 130 ±100 pC. The maximum achiev- 
able electron beam energy was found to increase from 70 
to 230 MeV for densities from 8 to 3 x 10^^ cm~^. Figure 
|3] shows the scaling of the electron energy with density. 
At low densities, electron beams were characterised by a 
high energy mono-energetic features, with a low energy 
dark current (fig. [2|l). For increasing densities, electron 
beams developed a more complex poly-energetic struc- 
ture and eventually a 100% energy spread. 
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FIG. 2: (a) Sample electron data. Electron charge per mrad 
and MeV (2D image), and radially (red) and spectrally (blue) 
integrated curves. The dashed lines indicate the detector limit 
and noise level, (b) Sample x-ray half-shadow data (2D im- 
age), experimental x-ray intensity integrated along the edge 
(black line) and modeled x-ray intensity using Fresnel diffrac- 
tion theory (dashed red, dotted green, dash-dotted blue lines). 

The propagation of the electrons is not affected by the 
magnet in the direction perpendicular to the dispersion 
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FIG. 3: Maximum electron energy (blue circles, left axis) and 
1/e^ electron bunch radius (red squares, right axis) as a func- 
tion of density. Horizontal error bars are ±10% due to un- 
certainty in density determination and vertical error bars are 
due to uncertainty in electron energy (±5%) and x-ray source 
size determination (95% confidence interval). 



of the spectrometer and thus the vertical electron beam 
divergence can be obtained from the ID dispersed signal 
(fig.|2^). Since beams are not necessarily mono-energetic, 
we integrate the dispersed spectrum along the energy axis 
to obtain an average RMS (root mean square) beam di- 
vergence. We find that the divergence shows little corre- 
lation with the density but may increase somewhat with 
higher densities. This was was observed before [rG", TT . 

It has been shown that electrons accelerated in a 
plasma wave produce copious amounts of synchrotron- 
like radiation, termed betatron x-rays [18] [19]. This is 
due to the electrons performing transverse oscillations 
whilst being accelerated along the laser axis in the fields 
of the plasma bubble, as illustrated in figure [l]3. The 
characteristics of the betatron x-rays are related to the 
properties of the electron beam and plasma wiggler. In 
particular the small x-ray source size can serve as a di- 
agnostic for the transverse electron bunch size. With 
our choice of laser parameters, the spot size is matched 
to the bubble radius wq = rij = 2^/a ■ c/up for densi- 
ties of 3.3 — 4.9 X 10^^ cm~^ and the pulse duration is 
matched to the bubble radius c ■ II = ^5 for densities of 
2.9 — 4.6 X 10^^ cm~^. If these conditions are fulfilled, the 
Wakefield accelerator operates in the bubble regime [4l. 
Dephasing then coincides with depletion and maximises 
the electron beam energy and thus the x-ray power at 
the end of the interaction. The source size of the x-rays 
therefore gives the transverse bunch size r± of the elec- 
tron beam at the end of the interaction. 

As shown in figure [l^, penumbral imaging is used to 
determine the betatron x-ray source size. A half-plane 
(thinly cleaved InSb crystal) is placed after the x-ray 
source, casting a half-shadow on the detector. The exper- 
imental intensity distribution integrated along the edge 
is shown in figure and can be modeled with Fresnel 
diffraction theory, yielding x-ray spectrum and source 



size [20]. Using a least-squares method, the RMS x-ray 
source radius (transverse electron bunch radius) is cal- 
culated with 95% confidence interval, as shown in figure 
|3j right axis. Much like the divergence, the x-ray source 
size also shows little correlation with density. 

Using equ.[l]and|2j the measured 7||, v±/v\\ and r± can 
be combined to obtain an upper limit for the RMS trans- 
verse normalized emittance, which is plotted in figure [4] 
Despite the fact that each measured quantity shows little 
correlation with density, the emittance gives a correlation 
of C = 0.84 with density. The associated p- value (prob- 
ability of getting a correlation as large as the observed 
value by random chance, when the true correlation is 
zero) is P < 10~^. A quadratic fit gives an R'^ = 0.83. 

Two effects will determine the scaling of the emittance. 
Firstly, for constant laser pulse duration, the emittance 
will grow at high densities, as the electron beam is more 
likely to interact with the laser pulse [l6l[T9]. This should 
become significant as c-t/, > r^, i.e. for densities > 2.9 — 
4.6 X 10^^ cm~^. This trend is reflected in our data (figure 
p|. Secondly, the emittance will scale with the bubble size 

— 1/2 

(X He and the transverse electron momentum spread ex 
ao [4 , leading to an increase at low densities. This trend 
is not reflected well by our data, which does not extend 
to densities much lower than where c ■ II — ^5 (shaded 
area in figure [4|. The fit to our data however suggests 
that the emittance is minimized at ng — 4 x lO^^cm"^, 
which is when wq = and c-tL = are matched. 

An analytic derivation for the bubble regime us- 
ing the FWHM quantities for {r±^pj_) gives the 
FWHM transverse normalized emittance as = 

0.033 aQ^^(ncrit/^e)^^^7 which computes to 0.87 
and 1.597rmmmrad for rig = 4 x 10^^ cm~^ and a = 
3.2 and 4.8 respectively [21]. The latter value is the 
self-consistent trapping threshold for the given density. 
The predicted emittance is in very good agreement with 
the measured RMS normalized transverse emittance o:^ 
0.5 TT mm mrad at ng = 4 x 10^^ cm~^ . 

In computing the measured emittance, we used the 
spectrally averaged divergence v_\_/v\\ and radially aver- 
aged mean energy 7|| . The radiated x-ray power per solid 
angle is ex 7^ which affects the measured x-ray source size 
and electron bunch size r±. This may motivate a weigh- 
ing of the 2D electron data with oc 7^ before computing 
averaged vj_/v\\ and 7||. When done, we find that both 
the absolute number and density trend of the emittance 
differs insignificantly from our figure [4| 

Comprehensive numerical modelling was performed to 
verify if such small emittances can occur in the bub- 
ble regime. Three-dimensional particle in cell simula- 
tions with the code OSIRIS [22] were carried out, in 
which a pulse with experimental parameters (ao = 3.0, 
t/, = 32 fs, (iFWHM = 10.8 /im) is focused 0.25 mm 
into the plasma. The electron plasma density profile in- 
creases linearly from zero to ng = 6 x 10^^ cm~^ in the 
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FIG. 4: Experimentally determined emittance as function of 
density (blue squares). Fit with a polynomial of 2nd order 
(red solid line), giving = 0.85 and 95% confidence interval 
(dashed red). The vertical grey lines indicate the density 
interval for which the spot size wq = rt and pulse duration 
c ■ tL = rb matches the bubble radius. 
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FIG. 5: (a) Electron trajectories obtained from the modeling, 
(b) Histogram of the betatron amplitudes for entire 3.7 mm 
jet (solid black) and for the last 1 mm (dashed red). 



first 0.5 mm, it is constant for 2.7 mm, and falls lin- 
early to zero in 0.5 mm. To save computational time, 
simulations were performed in a relativistically boosted 
frame (7 = 5) [23 . The simulation box corresponds to 
3.5 X 0.11 X 0.11 mm^ in the laboratory frame. A total 
of 9 X 10^ particles were pushed for 6 x 10^ iterations. 
The resolution in the laser propagation direction (in the 
boosted frame) is k^Az = 0.12 and kpAr± = 0.16 in the 
transverse directions. The modeling predicts a polyen- 
ergetic electron beam with peaks at 210 and 300 MeV, 
an average energy of 240 MeV and an almost symmetric 
profile with divergence v_\_/v\\ 4.6 mrad. As the code 
tracks the particle trajectories we find that electrons are 
injected throughout the interaction and all of them per- 
form oscillating trajectories as depicted in figure [5^. We 
analyze the oscillation radius for each electron and oscil- 
lation to obtain a histogram of oscillation radii, shown 
in figure ^jp. From the histogram we obtain an average 
oscillation amplitude or bunch radius of r_\_ = 0.9 and 
1.0 /im respectively, depending on whether we average 



over the entire length of the trajectory or the last 1 mm. 
This computes to a transverse normalized emittance of 
e'j_ = 0.6 — 0.7 TT mm mrad which is in good agreement 
with the measurement and analytical calculation. 

We have proposed and used a technique to measure the 
emittance of a LWFA electron beam. The method infers 
the transverse electron bunch size in the bubble from the 
measurement of the betatron x-ray source size. With this 
method we find that the transverse normalised emittance 
of a 230 MeV electron beam in the bubble regime can be 

< 0.5 TT mm mrad similar to conventional accelerators. 
The method can be applied to multi-GeV electron beams 
anticipated from LWFA in the near future, where stan- 
dard single shot emittance measurements have limited 
accuracy and eventually fail. Our method will be cru- 
cial to better diagnose LWFA and help to progress them 
towards applications. 

This work was supported by U.S. NRC grant no. 38- 
09-953. The work of J.L.M., R.F., J.V. and L.O.S. is 
partially supported by FCT (Portugal). The simulations 
were performed on the 1ST Cluster (Portugal) and on 
Jugene (Germany). 



[1] T. Tajima et al., Phys. Rev. Lett. 43, 267 (1979). 
[2] G. Z. Sun et al., Physics of Fluids 30, 526 (1987). 
[3] A. Pukhov et al.. Plasma Phys. Control. Fusion 46, B179 
(2004). 

[4] W. Lu et al., Phys. Rev. Spec. Top.- Accel. Beams 10, 

061301 (2007). 
[5] S. P. D. Mangles et al.. Nature 431, 535 (2004). 
[6] J. Faure et al.. Nature 431, 541 (2004). 
[7] C. G. R. Geddes et al.. Nature 431, 538 (2004). 
[8] L. E. Collins et al., Nucl. lustrum. Methods 26, 157 

(1964). 

[9] Y. Yamazaki et al., Nucl. lustrum. Methods Phys. Res. 
Sect. A-Accel. Spectrom. Dect. Assoc. Equip. 322, 139 
(1992). 

[10] S. Fritzler et al., Phys. Rev. Lett. 92, 165006 (2004). 
[11] E. Brunetti et al., Phys. Rev. Lett. 105, 215007 (2010). 
[12] C. M. S. Sears et al.. Physical Review Special Topics - 

Accelerators and Beams 13, 092803 (2010). 
[13] W. P. Leemans et al., Nat. Phys. 2, 696 (2006). 
[14] S. Kneip et al., Phys. Rev. Lett. 103, 035002 (2009). 
[15] K. T. McDonald et al., in Frontiers of Particle 

Beams: Observation, Diagnosis and Correction, edited 

by M. Month et al. (Springer, 1989), p. 122. 
[16] S. P. D. Mangles et al., Phys. Rev. Lett. 96, 215001 

(2006). 

[17] S. Kneip et al.. Plasma Phys. Control. Fusion 53, 014008 
(2011). 

[18] A. Rousse et al., Phys. Rev. Lett. 93, 135005 (2004). 
[19] S. Kneip et al., Phys. Rev. Lett. 100, 105006 (2008). 
[20] S. Kneip et al., Nat. Phys. 6, 980 (2010). 
[21] A. G. R. Thomas, Phys. Plasmas 17, 056798 (2010). 
[22] R. A. Fonseca et al. (Springer, Heidelberg, 2002), vol. 
2331, p. 342. 

[23] S. F. Martins et al., Nat. Phys. 6, 311 (2010). 



